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Abstract 

We analyze heating and cooling processes in accretion disks in bina- 

• 

ries. For realistic parameters of the accretion disks in close binaries (M ~ 
10~^^ 10~'^ Mq/ year and a ~ 10^^ 10~^), the gas temperature in the 
outer parts of the disk is ~ 10^ K to 10^ K. Our previous gas-dynamical 
studies of mass transfer in close binaries indicate that, for hot disks (with 
temperatures for the outer parts of the disk of several hundred thousand K), 
the interaction between the stream from the inner Lagrange point and the 
disk is shockless. To study the morphology of the interaction between the 
stream and a cool accretion disk, we carried out three-dimensional modeling 
of the flow structure in a binary for the case when the gas temperature in the 
outer parts of the forming disk does not exceed 13 600 K. The flow pattern 
indicates that the interaction is again shockless. The computations provide 
evidence that, as is the case for hot disks, the zone of enhanced energy re- 
lease (the "hot line") is located beyond the disk, and originates due to the 
interaction between the circum-disk halo and the stream. 



1 Introduction 

In 1999-2002, we developed a three-dimensional, gas-dynamical model and used 
it to study the flow patterns in binary systems ifTUTill. These studies indicate that 
the flow structure is substantially affected by rarefied gas of the intercomponent 
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envelope. In particular, a self-consistent solution does not include a shock in- 
teraction between the stream from the inner Lagrange point Li and the forming 
accretion disk (a "hot spot"). The region of enhanced energy release (the "hot 
line") is located beyond the disk and is due to the interaction between the enve- 
lope and the stream. However, these solutions were obtained for temperatures of 
the outer parts of the accretion disk of 200 000 - 500 000 K. To check if this 
behavior is universal, the morphology of the flow must be considered for various 
disk temperatures. 

First, we will study here the interval of plausible temperatures of accretion 
disks in close binaries. In Section [2l based on an analysis of heating and cooling 
in accretion disks, we will show that, for realistic parameters of the disks in close 

binaries, (M~ 10"^^ -r- IO'^Mq /year and a ~ 10^^ ^10^^), the gas temperature 
in the outer parts of the disk is between 13 600 K and ~ 10^ K. This implies that 
cool accretion disks can form in some close binaries. 

Second, we will consider the morphology of the interaction between streams 
of matter and cool accretion disks in semi-detached binary systems (Sections [3] 
andlH). The basic problem here is whether the interaction between the stream and 
the disk remains shockless, as was shown for relatively hot disks [fn[3ll4l [811141. 
Section [5] presents our main conclusions and a physical basis for the universal 
nature of the shockless interaction between the stream and disk. 

2 Heating and Cooling in Accretion Disks 

In this Section, we consider the temperature of an accretion disk for various ac- 
cretion rates, i.e., the dependence T(M). 

2.1 Basic Equations 

The vertical structure of an accretion disk is specified by the balance between the 
vertical component of the gravitational force and the (vertical) pressure gradient, 
which, in turn, is specified by the balance between heating and cooling of the 
gas. The heating is associated with viscous dissipation of kinetic energy, and also 
with bulk radiative heating, which, in turn, is specified by the radiation of the 
central object. Cooling is brought about by several mechanisms: bulk radiative 
cooling, radiative heat conduction, and convection. Assuming that advective terms 
and terms associated with adiabatic heating or cooling are small, the steady-state 
energy equation 

g+-g- = 

can be written as follows. 
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(1) For the optically thin case, when is specified by bulk radiative heating 
and viscous heating and is determined by bulk radiative cooling, 

Q+,,(p, T)+n^. (r(T, T^,) - A(T)) = . (1) 

Here, T{T,T^ii) is the radiative-heating function, which depends on the gas 
temperature T and the temperature of the central object Tyjd, A(T) is the radiative- 
cooling function, and Q^isdp, T) is the viscous heating. 

(2) For the optically thick case, is specified by viscous heating, while Q~ 
is specified by radiative heat conductiorQ and convection in the vertical direction: 

QvisciP,T) ~ -Q^=0. (2) 

Here, Frad and Fconv are the radiative and convective energy fluxes. To determine 
the functions in ([T]) and 0, we will need 

- the equation of continuity 

— M = 2tt / r ■ p ■ Vr dz = const , (2) 



- the equation of angular-momentum balance A = r'^Qx in the radial direc- 
tion: 

9-,(r-P--r-X) = g-^[u-p-r — j , (3) 
from which it follows that 

\vr\ = -ly ■ n'j^ ■ ■ r-^ ^ ly ■ r-\ (4) 



- and the viscous heating 



Here, Mis the accretion rate, = yGM/r^ the angular velocity of the Ke- 
plerian rotation of the disk, M the mass of the central object, G the gravitational 
constant, p - the density, Vr the radial velocity, and the u - coefficient of kine- 
matic viscosity. Note that the molecular viscosity cannot provide the necessary 

'We neglect molecular heat conduction since it is very small compared with the radiative heat 
conduction. 
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dissipation, and dissipation processes are usually considered to be associated with 
turbulent or magnetic viscosity. 

To determine the vertical pressure gradient, we will use the equation of hydro- 
static balance in the vertical direction 



1 dP Of GM 



Vl\z, (6) 



p dz dz \ ^Jr"^ -\- z"^ J 

as well as the equation of state of an ideal gas with radiation 

P = pTZT + VsaT^ . 

Here, P is the pressure, T the temperature, IZ the gas constant, and a the radiation 
constant. All equations are given in cylindrical coordinates, (r, z). 



2.2 The Solution Method 

To determine the dependence T{M), we will use ([2l) and (HI) together with the 
expression for the viscosity coefficient . We will use the formula for suggested by 



Shakura [fTSlI , u = aCgH, where H is the height of the disk and ~ y' TZT + Y ^oT^ j p 
is the sound speed. If we neglect the z dependence of the density and use p av- 
eraged over the height (further, we will denote this quantity simply as p), the 
integration of Q yields the height of the disk H: 

We will determine Cg from the temperature in the equatorial plane of the disk, 
z = 0. This approach is sufficiently correct for our purposes due to the uncertainty 

in the parameter a. As a result, we obtain an equation relating M, T|^^q, and p 
for the specified r and a, 

M=2TT-a-QK^-p-cl = 2n-a-QK^- (iZTp^/^ + y^aT^p-^/^f^^ . (7) 

This equation reduces to a cubic equation in the variable p^^^, and its solution has 
two branches: one with a negative real root and two complex ones, and one with 
three real roots, one of which is negative. Only positive real roots for the density 
are physically meaningful. For such roots to exist, the following condition must 
be satisfied: 

• VSn -Vn-a-a- T^''^ 

M> -2 . (8) 

which yields the minimum accretion rate for the given T, r, and a. When deriving 
this condition, we used the equation of state taking into account the radiation 
pressure. 
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This estimate can also be written in the form 



/ . X 2/9 



where R^d = lO^cm is the radius of the accretor (white dwarf). 

Let us consider the condition ([8]) for the outer parts of the accretion disk. Let 
us take r = A/ 5, where A is the distance between the components of the binary 
A = 1A2Rq), which corresponds to situation for the dwarf nova IP Peg; as a 
result, we obtain 

• / T \^/^ / a \ 

^>10 T7^ 77T M^/year. (9) 



10° K/ VO.l 

If © is satisfied, the roots of Eq. ([7]) relating p, T, and M for a given r and a can 
be written 

p = (7^T)-3/4 I ^\ 3i^3 / i/3arcsin { v^aT^/^^""" 



(for simplicity, we have omitted numerical factors — 1). The first of these 
corresponds to disks with dominant radiation pressure (/3 = ^/^aT^ / pTZT > 1) 
and the second to disks with dominant gas pressure {(3 < 1). 

These formulas describe the two branches of the two-parameter dependence 

p{M, T). To calculate the dependence T{M), we must use the additional heat bal- 
ance equations ([T])-©. As follows from Section [2?T1 the form of ©-© depends 
on the optical depth of the disk, which, accordingly, must be calculated. 



2.3 Optical Depth 

The optical depth r is specified by the product of the absorption coefficient k, the 
density, and the geometrical depth of the layer I: t = k ■ p ■ I. For disk accretion, 
the basic parameter is the ratio of the geometrical depth of the layer where r = 1 
and the height of the disk: /H. After simple manipulation, we obtain 



10 ' 10 ' 10 ° 
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Figure 1: The k(T) dependence for n = lO^^cm '^,n = 10^ ''cm = lO^^cm ^, 
n = W^^cm^^, and n = W^^cm^^ (top to bottom) ifTSl . 

The absorption coefficient k displays a complicated dependence on T and 
p (and also on the degree of ionization, chemical composition, etc.). Here, we 



adopted the simple approximation for / 


<T,p)|[l6 
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«4 , 


K4 = 


= 0.348 . 





According to [18J, these four subregions correspond to scattering on molecular 
hydrogen, scattering on atomic hydrogen, free-free and free-bound transitions, 
and Thompson scattering. The boundaries of the sub-regions, i.e. the transitions 
from one expression to another, are specified by the equality of the k values cal- 
culated from these expressions. Figure [T] presents the dependences of k on T and 
p. We can see regions with dn/dT > 0, where thermal instability can develop 
when the dependence between the surface density and the disk temperature forms 
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Figure 2: Solution of (|7]) for disks with dominant gas pressure in the (T, M) plane for 
a = 0.1 and r = A/5. The horizontal thick shading indicates optically thick disks, and 
vertical shading optically thin disks; the solid line marks the border between these regions. 
The dashed line corresponds to condition Q for the existence of the solution of there 
is no solution below this line. 

an S curve in the (S, T^ff) plane. Thermal instability is often invoked to explain 
the phenomenon of dwarf novae (see, for example, [UQHIQID : however, it is clear 
that this can occur only for sufficiently cool disks. 

Let us return to (|7]), taking a = 0.1 and r = A/h and considering disks with 
dominant gas pressure, for which (3 = l/3aT^/ pTZT < 1. The shaded region in 

the (T, M) plane in Fig.[2]corresponds to all possible solutions for these disks. The 
dashed line corresponds to condition Q for the existence of a solution for (|7]) - 
there is no solution below this line. The solid line indicates the boundary between 
the optically thick and optically thin solutions: the horizontal shading marks the 
region of optically thick disks, while the vertical shading marks optically thin 
disks. Figure [3] presents a similar pattern for disks with dominant radiation pres- 
sure {/3 > 1). 

We can see from Fig.[2]that, for realistic values M G [10^^^, 10^'^]MQ/year, 
disks with dominant gas pressure are mainly optically thick, though solutions cor- 



7 



10 " 10 ' 10 

T (K) 



Figure 3: The same as Fig.|2]for disks with dominant radiation pressure. 

responding to optically thin cool disks are possible for small M. It follows from 
Fig. [3] that disks with dominant radiation pressure are mainly optically thin; opti- 

cally thick hot disks can exist only for high M. 

2.4 Optically Thick Disks 

In Section [Z2l we derived Eq. (|7]), which relates M, T, and p for given a, r and 
a . Using the supplementary heat-balance equations ([T])-©, we can reduce the 

number of unknowns and obtain the desired relation between M and T. 

The vertical temperature distributions in optically thick disks are described by 
the equation of radiative heat conduction with a source due to viscous heating Q, 
which can be written in the form 




where e is the specific internal energy and c is the velocity of light. To solve (fTTI) . 
we must specify boundary conditions. Due to the symmetry of the problem, the 
temperature derivative in the equatorial plane must be zero; i.e., T'|^^q = 0. 
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The temperature at the upper boundary of the disk is specified by the condition 
r(T,,T^d) = A(T,). Though the functions T{T,T^d) and A(T) are complex, 
they are known and can be found in the literature (see, for example, [|2T1 - [231 ). The 
temperature derived by equating these functions (for a temperature of the central 
object (white dwarf) of T^^ = 70 000 K) is T{H) = = 13 600 K. 

The solution of (fTTI) enters a stationary regime when the characteristic heat- 
conduction time 

is comparable to the time for viscous heating 

Note that (fTTI) can be integrated analytically in the steady-state case. Let us denote 
U = T^, = T^, Uq = U\^^Q and again assume that does not depend on z. 
Then, 

d I 1 d (ac. 



After integrating over z, we obtain 

1 d fac \ 2o 
npdz \ 6 J 

The integration constant is equal to zero, since IJ'\^^^ = 0. For convenience, we 
will transform this last equation to the form 

IdU dB 3 2 2r. 

-— = — = p ac^ilKZ, 

K dz dz ac 

, , dB 
where the function BiJJ) is determined from the differential equation = 

and can be written in an analytical form if p is fixed. Integrating this 



<U,p) 

last equation over z, we obtain 



or, for z = 



B{Uo) = B{U,) + ^p^acl^KH^ ■ 
2ac 



Using the expressions 



2 I'DtA/^: I 1/ 

Cs = l^Uo' + 73 — 
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Figure 4: Solution of (ITTI) for an optically thick disk for a = 1 and r = A/b (asterisks). 
Solutions of ([T2l ) taking into account convection are labeled by squares. The dashed line 
is a lower bound for the region in which there exist solutions of d?), and the solid line 
separates the regions of optically thin and optically thick disks. 

we obtain the algebraic equation for Uq 

This equation implicitly specifies the dependence Uo{p); i.e., T(p). Expressing M 

in terms of p and T, we can derive the dependence M(p) = MiT{p), p), which 

yields the dependence T[M) in parametric form. Formally, the resulting solution 
can also exist in optically thin regions; however, given the adopted assumptions, 
these points can be rejected. 

Let us consider a graphical representation of the solution derived. Figure |4] 

presents the dependence T(M) for a = 1 and r = A/5, marked by asterisks. 
As in Fig. [2]-[3l the dashed lines bound from below the domain in which there 
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Figure 5: The same as Fig. |4] for a = 0.1. 



exist solutions of ([7]), while the solid line separates the domains of optically thin 
and optically thick disks. Figures [5]-[6] display the solutions for a = 0.1 and 
a = 0.01, respectively. Figure [6] presents the accretion rate as a function of the 
disk thickness. We can see that all the obtained disks are geometrically thin; i.e., 
i7 < r. 

Radiative heat conduction is not the only mechanism for heat transfer into 
optically thin regions. Under certain conditions, convection can also play a sub- 
stantial role. Neglecting the radiation pressure, the convective flux can be written 
in the form [I24ll25]l 



Cp ■ p ■ 



AVT 




1/2 

T 

Cp 



4 

dS_ 

dz 



(AVT) 



3/2 



Here, cp is the heat capacity at constant pressure, S = TZ ■ ln(T^/^/ p) the spe- 
cific entropy, g = —Vl]^z the gravitational acceleration, and / the mixing length, 
taken to be / = aH. To determine the vertical temperature distribution taking 
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Figure 6: The same as Fig.|4]for a = 0.01. 
convection into account, we must solve the equation 

with the same boundary conditions as for (fTTI) . Equation (fT2l) does not admit 
a simple analytical solution, and we solved this equation numerically using the 
method of establishment. The solution is denoted by the squares in Figs|4]-[6l We 
can see that convection plays a significant role only when a ~ 1. 

Summarizing, we can assert that, in the optically thick case with small M, the 
disk displays the constant temperature T = = 13600° K, while the temperature 

increases as T oc M^^'^ at larger values of M . Thus, for realistic parameters 

of the accretion disks in close binaries, M 10~^^ 4- 10~'^MQ/year and a ~ 
10"^ ^ 10"^, the gas temperature in the outer parts of the disk (r ~ A/5 A/10) 
is 10^ K to ~ 10^ K. 

Solving (fTTI) for various r, we can also calculate the dependences T(r) and 
p{r). The calculations indicate that T oc r^°-^ and p oc r~^-^, which is consistent 
with the dependence T oc r^'^/^ obtained by Shakura and Sunyaev [26J. 
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Figure 7: A disk witli dominant gas pressure. The solid lines indicate possible states 

of the disk in the H/r — M plane for r = A/b. Solutions of (ITTI) taking into account 
radiative heat conduction and viscous heating are shown by the lines with asterisks for 
a = 1, a = 0.1, a = 10^^, a = 10^^ (top to the bottom). Solutions of ([12]) taking into 
account radiative heat conduction, convection, and viscous heating are shown by the lines 
with squares for a = 1. The dashed lines bound from below regions in which the solution 
of © can exist for a = 1, a = 0.1, a = 10~^, a = 10^^ (top to bottom). 

2.5 Optically Thin Disks 

In this case, the temperature of the disk is specified by the balance between ra- 
diative heating r(T, Tyj^) and viscous heating ([5]), on the one hand, and radiative 
cooling A(T), on the other. The heat-balance equation ([T]) can be written 

pac^^^K + p'-ml- (r(T, T^,) - A(T)) = , 

which can be reduced to the quadratic equations in p 

a ■ {pUT + VaaT^) ■ + p" ■ mf ■ (r(T, T^,) - A(T)) = . 

The solution of this equation for specified r and a yields the dependence p(T), 

and thus T(M). Formally, this solution can also exist in optically thick regions, 
however these points were rejected by virtue of the adopted assumptions. 
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Figure 8: Solutions for an optically thin disk for a = 1, a = 0.1, a = 10^^, a = 10^^ 
(top to bottom) and r = A/20 (asterisks). The dashed lines bound from below the domain 
in which there exists a solution of Q; the solid lines separate the regions for optically thin 
and optically thick disks. 

It is shown in Section [23] that disks in which gas pressure dominates are pri- 
marily optically thick, and solutions that correspond to optically thin disks are 

possible only for small M. Disks in which radiation pressure dominates are pri- 
marily optically thin. The domination of radiation pressure is possible only in the 
inner parts of the disk; therefore, we will adopt r = A/20 for the further analysis. 
For the typical dwarf nova IP Peg, this corresponds to five radii of the accretor 
(white dwarf). 

Figure [8]presents the results of our calculations; the asterisks denote the T(M) 
dependences for a = 1, a = 0.1, a = 10~^, a = 10~^ (top to bottom), and 
r = A/20. The disks obtained in these solutions are geometrically thick, H c:^ r. 
Note that the initial assumptions of the model restrict its applicability: it is suitable 
only for geometrically thin disks, and the solutions for geometrically thick disks 
are purely formal. 
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3 The Model 



We described the flow structure in a binary system using a system of gravitational 
gas-dynamical equations taking into account radiative heating and cooling of the 
gas for the optically thin case: 



Here, as usual, p is the density, v = (u, v, w) the velocity vector, P the pressure, 
e the internal energy, $ the Roche gravitational potential, nip the proton mass, 
and r(T, Tu,d) and A(T) the radiative heating and cooling functions, respectively. 
The system of gas-dynamical equations was closed with the Clapeyron equation 
P = ('J — l)pe, where 7 is the adiabatic index. We took the parameter 7 to be 5/3. 

Our main goal here is to study the morphology of the interaction between the 
stream and the cool accretion disk. It follows from Section |2] that the outer parts 

of the accretion disk can be cool for small M and, in particular, in the case of an 
optically thin disk. The system of equations (fT3l) enables us to carry out three- 
dimensional modeling of the flow structure in a binary within our formulation of 
the problem. In the model, the temperature of the disk is 13 600 K. 

We solved this system of equations using the Roe-Osher method [[T4ll27ll28l . 
adapted for multiprocessing computations via spatial decomposition the compu- 
tation grid (i.e., partitioning into subregions, with synchronization of the bound- 
ary conditions) [|29ll . We considered a semi-detached binary system contain- 
ing a donor with mass M2 filling Roche lobe and an accretor with mass Mi. 
The system parameters were specified to be those of the dwarf nova IP Peg: 
Ml = I.O2M0, M2 = O.5M0, A = 1.42i?0. 

The modeling was carried out in a non-inertial reference frame rotating with 
the binary, in Cartesian coordinates in a rectangular three-dimensional grid. Since 
the problem is symmetrical about the equatorial plane, only half the space occu- 
pied by the disk was modeled. To join the solutions, we specified a corresponding 
boundary condition at the lower boundary of the computation domain. The accre- 
tor had the form of a sphere with radius 10^"^ A. All matter that ended up within 
any of the cells forming the accretor was taken to fall onto the star. A free bound- 
ary condition was specified at the outer boundaries of the disk - the density was 



< 




(13) 



= -pv grad $ + p2m;2 (r(T, T^, 



:)-A(T)). 
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constant (p^ = lO^SpLi)^ where is the density at the point Li, the tempera- 
ture was 13 600 K, and the velocity was equal to zero. The stream was specified 
in the form of a boundary condition: matter with temperature 5800 K, density 
PLi = 1-6 X 10^^(7/cm^ and velocity along the x axis = 6.3km/s was injected 
into a zone around Li with radius COMA. For this rate of matter input into the 
system, the model accretion rate was ~ 10^^"^ Mq/ year. 

The size of the computation domain, 1.12A x 1.14A x 0.17A, was selected 
so that it entirely contains both the disk and stream, including the point Li. The 
computation grid with 121 x 121 x 62 cells was distributed between 81 processors, 
which constituted a two-dimensional 9x9 matrix. 

To increase the accuracy of the solution, the grid was made denser in the zone 
of interaction between the stream and disk, making it possible to resolve well the 
formed shock wave. The grid was also denser towards the equatorial plane, so that 
the vertical structure was resolved, even for such a cool disk. 

We used the solution obtained for a model without cooling as the initial con- 
ditions [fT2ll . The model with cooling was computed during approximately five 
revolutions of the system, until the solution became established. The total time of 
the computations was ~ 1000 h on the MBCIOOOA computer of the Joint Super- 
computer Center (JSC). 

4 Computation Results 

Figures l9lto [T3]present the morphology of gas flows in the binary. Figure |9] shows 
the density and velocity vector distributions in the equatorial plane of the system 
(the XY plane), while Figs. [10] and [TT] present density contours in the frontal 
(XZ) plane and in the YZ plane containing the accretor and perpendicular to the 
line connecting the binary components. In spite of the small height of the forming 
accretion disk, use of the JSC parallel-processing computers made it possible to 
resolve its vertical structure (the outer parts of the disk were covered by 15 grid 
cells, and the inner parts by no fewer than 3 cells). Figure [I2l gives an enlarged 
view of the density and velocity vector distributions in the zone of interaction 
between the stream and the outer edge of the disk (the area in the shaded rectangle 
in Fig.|9l). Figure [13] presents the so-called texture a visualization of the velocity 
field in the zone of interaction between the stream and disk, constructed using the 
Line Integral Convolution procedure (LIC) |[30l . 

According to our considerations in |[8l[T4l, the gasdynamical flow pattern in 
a semi-detached binary is formed by the stream from Li, the disk, a circumdisk 
halo, and the intercomponent envelope. This subdivision is based on physical dif- 
ferences between these elements of the flow structure: (1) if the motion of the gas 
is not determined by the gravitational field of the accretor, it forms the intercom- 
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Figure 9: Contours of constant density and velocity vectors in the equatorial plane XY 
of the system. The shaded rectangle indicates the zone of interaction between the stream 
and disk, shown in Figs[T2]and[T3l The point Li and the direction towards L3 are marked. 
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Figure 10: Density contours in the frontal plane XZ of the system. 
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Figure 1 1 : Density contours in the plane YZ containing the accretor and perpendicular 
to the line connecting the binary components. 

ponent envelope; (2) if the gas makes one revolution around the accretor, but later 
mixes with the initial stream, this gas does not become part of the disk, instead 
forming the circum-disk halo; (3) the disk is formed by that part of the stream 
that loses its momentum and moves towards the center of gravity after entering 
the gravitation field of the accretor, rather than interacting with the stream. 

In this framework, let us consider the morphology of the flow when the tem- 
perature decreases to 13 600 K over the entire computation domain due to cooling. 
Figure |9] indicates that, in this case, the intercomponent envelope is formed pri- 
marily in the vicinity of L3, and does not affect the solution substantially. We can 
see from Figs. I9l and [TO] - [TT] that the circum-disk halo is pressed against the disk, 
and its density increases sharply towards the disk edge. 

Figures [121 and [13] show that, in the cool-disk case, the interaction between the 
circum-disk halo and the stream displays all features typical of an oblique collision 
of two streams. We can clearly see two shock waves and a tangential discontinuity 
between them. The gases forming the halo and stream pass through the shocks 
corresponding to their flows, mix, and move along the tangential discontinuity 
between the two shocks. Further, this material forms the disk itself, the halo, and 
the envelope. 

The solution for the cool case displays the same qualitative characteristics as 
the solution for the case when the outer parts of the disk are hot: the interaction 
between the stream and disk is shockless, a region of enhanced energy release is 
formed due to the interaction between the circum-disk halo and the stream and is 
located beyond the disk, and the resulting shock is fairly extended, which is par- 
ticularly important for explaining the observations. However, unlike the solution 
with a high temperature in the outer regions of the disk [|Tlll2l[14l, in the cool 
case, the shape of the zone of shock interaction between the stream and halo is 
more complex than a simple "hot line". This is due to the sharp increase of the 
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Figure 12: Contours of constant density and velocity vectors in the zone of interaction 
between the stream and disk (the shaded rectangle in Fig.|9ll. 



19 



-0.14 




Figure 13: Visualization of the velocity field in the zone of interaction between the stream 
and disk (the shaded rectangle in Fig.|9l). 
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halo density as the disk is approached. Those parts of the halo that are far from 

the disk have low density, and the shock due to their interaction with the stream 
is situated along the edge of the stream. As the halo density increases, the shock 
bends, and eventually stretches along the edge of the disk. 

5 Conclusions 

Our analysis of the basic processes of heating and cooling in accretion disks in 
binaries has shown that, for realistic parameters of the accretion disks in close 

binary systems (M ~ 10~^^ 10~'^MQ/year and a ~ 10~^ 10~^), the gas 
temperature in the outer parts of the disk is 10"^ K to ~ 10^ K. 

Previously, we carried out three-dimensional simulations of the flow structure 
in close binaries for the case when the temperature of the outer parts of the accre- 
tion disk was 200-500 thousand K. Those solutions showed that the interaction 
between the stream from the inner Lagrange point and the disk was shockless. To 
determine the generality of the solution, the morphology of the flow for different 
disk temperatures must be considered.We have presented here the results of sim- 
ulations for the case when cooling decreases the temperature to 13 600 K over the 
entire computation domain. 

Our analysis of the flow pattern for the cool outer parts of the disk confirms 
that the interaction between the stream and disk is again shockless. The computa- 
tions indicate that the solution for the cool disk case displays the same qualitative 
features as in the case when the outer parts of the disk are hot: the interaction 
between the stream and disk is shockless, a region of enhanced energy release 
formed by the interaction between the circum-disk halo and the stream is located 
beyond the disk, and the shock wave that is formed is fairly extended, and can be 
considered a "hot line" . The cool solution demonstrates the universal character 
of our previous conclusions that the interaction between the stream and disk in 
semidetached binaries is shockless. 
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